Reconstruction of field theory from excitation spectra of defects 
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We show how to reconstruct a field theory from the spectrum of bound states on a topological 
defect. We apply our recipe to the case of kinks in 1+1 dimensions with one or two bound states. 
Our recipe successfully yields the sine-Gordon and A<^> 4 field theories when suitable bound state 
spectra are assumed. The recipe can also be used to globally reconstruct the inflaton potential of 
inflationary cosmology if the inflaton produces a topological defect. We discuss how defects can 
provide "smoking gun" evidence for a class of inflationary models. 

PACS numbers: 03.65 



I. INTRODUCTION 

Often one is interested in diagonalizing a known oper- 
ator. This is an eigenvalue problem. The reverse prob- 
lem, where the eigenvalues are known but the operator 
is not, is also of interest. This is the "inverse scatter- 
ing problem" 0, 0] (also see e.g. Jjj). Indeed it is the 
basis of the famous question: "Can one hear the shape 
of a drum?" 0,0- In other words, can the frequencies 
(eigenvalues) of the sound emitted from the drum be used 
to determine the shape (operator) of the drum? In the 
quantum mechanics context, the equivalent problem is of 
reconstructing the Schrodinger potential, U(x), from the 
energy spectrum. Such inverse scattering problems are 
of interest in a wide variety of applications. 

We will be interested in a class of problems where U 
itself originates from some field interactions. Our aim is 
not just to reconstruct the Schrodinger potential, U, but 
the underlying field theory interactions that led to that 
particular Schrodinger potential. The eigenvalue data 
that we will use to reconstruct the field theory is the ex- 
citation spectra of any topological defects that might be 
present in the system. The method is widely applicable 
since topological defects are present in a large class of 
systems. They are routinely observed in condensed mat- 
ter systems, and are inevitable in high energy particle 
physics models. Topological defects relevant for high en- 
ergy particle physics are also relevant for the very early 
universe. Hence the excitation spectra of these topologi- 
cal defects can be a window to the early universe, and in 
particular, to inflationary cosmology 1 . 

A way to understand the usefulness of topological de- 
fects in a reconstruction effort is that the core of a topo- 
logical defect contains a different phase of the theory and 
so the core holds non-perturbative information about the 
field theory. Therefore the excitation spectrum of the de- 



fect leads to a global reconstruction of the field theory. 
This is in contrast to a local reconstruction that is possi- 
ble by perturbative methods such as particle scattering 2 . 

In what follows, we assume that we are given a topo- 
logical defect. In fact, for actual calculations, we will ex- 
clusively work with kinks in 1+1 dimensions since this 
example illustrates the central ideas and also because 
there are extensive techniques in 1+1 dimensions that are 
not available in higher dimensions. (Although spherically 
symmetric systems in higher dimensions can be reduced 
to the 1+1 dimensional case.) 

In the next section we will outline the recipe for the re- 
construction and then in Sec. IIIII explicitly work through 
two specific examples with kinks in 1+1 dimensions. Our 
scheme yields the sine-Gordon field theory when there is 
only one bound state (the translation mode) on the kink. 
When the kink has two bound states, we give implicit ex- 
pressions for the field theory. For a specific choice of the 
eigenvalues of the two bound states, we obtain the A</> 4 
field theory. In Sec. I1VI we discuss the reconstruction 
problem in the inflationary context. Then we qualita- 
tively discuss the bound state spectrum on kinks made 
from the inflaton field. In particular, we look for signa- 
tures of the flatness of the inflaton potential in the spec- 
trum. Readers only interested in inflation can proceed 
directly to Sec. |IV| since the discussion there is largely 
independent of Sec. ITU and IIIII 



II. GENERAL RECONSTRUCTION SCHEME 



We assume that we know the spectrum of energy eigen- 
values, {w 2 }, for the bound states on a topological de- 
fect. For the time being, we will only focus on the bound 
states. The price we pay is that there is then a huge de- 



1 Even if the abundance of topological defects is suppressed in the 
cosmos, the inevitability of defects in particle physics implies 
that they can (in principle) be produced in accelerators. 



2 In the inflationary case, astrophysical data can be used with 
certain assumptions to reconstruct the inflaton potential. This 
reconstruction too is local since astrophysical observations only 
probe a very limited range of relevant scales. 
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generacy in the reconstruction. By using scattering infor- 
mation and further physical input it may be possible to 
eliminate the degeneracy but we shall not be discussing 
this issue here. 

The general reconstruction scheme is as follows. 

(i) Prior to an analysis, a theoretical framework is 
needed and the appropriate field content should be 
known. In this discussion we will only consider a 
single scalar field (f> in one dimension, with standard 
form of the Lagrangian: 



L = 2^ 



- V{cf) 



(1) 



(ii) The Schrodinger equation that determines the ex- 



citation spectrum is: 

U(x) 



dx 2 



4>r, 



where 



U(x) = V"{fo{x)) 



(2) 



(3) 



and 4>o(x) is the (unknown) profile function of the 
defect - the classical defect solution. We would like 
to determine the potential U(x). Inverse scattering 
methods have been developed precisely to solve this 
problem. The answer, however, is not unique, es- 
pecially if only bound state spectra are taken into 
account. However, additional theoretical input can 
possibly reduce degeneracies. For example, if some 
interactions are known by perturbative methods, 
the information might be useful to break some of 
the degeneracy 3 . 

(hi) Once we have U(x), we find the "translation mode" 
by solving the zero eigenvalue Schrodinger problem: 



-dx-* +U{x) 



flt = 



(4) 



(iv) The translation mode is simply related to the defect 
profile functions by differentiation: 



dx 



(5) 



Hence, we can integrate the translation mode to 
obtain </>o(x). 

(v) Next we invert <f>o(x) to obtain x(<fio). 



3 U(x) may also be reconstructed using scattering data i.e. not 
just the eigenvalues but the scattering amplitudes. For example, 
in the Born approximation, the scattering amplitude is directly 
related to the Fourier transform of U(x). 



(vi) The equation of motion for the defect is: 



£ + V(A,) = o 



dx 2 



(6) 



where prime refers to differentiation with respect 
to the argument. Combining with Eq. JSJ) we see 
that: 



dipt 
dx 



(7) 
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(vii) Finally, an integration yields the desired symmetry 
breaking potential: 



V{<h) = / d<p v'(M 



(8) 



In the 1+1 dimensional case described above, there is 
actually a slight shortcut that is available to us. In this 
case, the Bogomolnyi equation holds: 



dx 



and so: 



(9) 



(10) 



x{4>o) 



There are two parts to this recipe. The first is the 
derivation of U (x) from the eigenvalue spectrum and the 
second is the derivation of V{<j>) from U{x). The part of 
the recipe starting with U(x) and constructing V((j>) as 
in Eq. (|10fl has also been used in earlier work 

The hardest step in this scheme is step (ii), the re- 
construction of U{x) from {w^}. In 1+1 dimensions (or 
for S-wave states in 3+1 dimensions), a simple general 
scheme to reproduce the bound state spectrum is given 
in Ref. @,H,11 The scheme employs the idea of super- 
symmetric quantum mechanics where the Hamiltonian 
operator can be factored, and yields a reflectionless po- 
tential with the desired bound state spectrum. We now 
summarize the scheme; a sketch of how the scheme is de- 
rived is given in Appendix [S] While it may seem that 
the reconstruction in terms of supersymmetric potentials 
is overly restrictive, we show in Appendix[5]that kinks in 
one spatial dimension always lead to a supersymmetric 
form of the potential U(x). We also show that this is 
can be true even when we have multi-component fields 
in more than one dimension. Thus the reconstruction 
scheme of Q is perfectly suited to our context. 

Suppose the bound state eigenvalues are labeled in de- 
scending order: J 2 > uif +1 and i = 1, N — 1 and the 
zero energy level is chosen so that the ground state eigen- 
value, Wat, is zero. Then a potential containing n of the 
highest bound states is: 



U n (x) = fl + f> n+ ul 



(11) 
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where the function f n (x) satisfies: 

f'n - fn + Un-l = U) n 



(12) 



III. EXAMPLES 



One bound state 



in terms of U n -\, the potential containing n — 1 of the 
highest bound states. If we write f n (x) = —w' n /w n , the 
equation for w n is: 



W„ + U n -iW n = UJ n W n 



(13) 



This equation will have two linearly independent solu- 
tions. If we further require that U n be even under parity: 
U n (-x) = +U n (x), then we need f n (—x) = -f n (+x), 
and w n {— x) = w n (+x). This condition eliminates one of 
the linearly independent solutions 4 . 

The solution to the n th Schrodinger equation: 



tp n = UJ n 1p„ 



is simply: 



1pn = 

Wn 



(14) 



(15) 



where a n is a normalization constant to be determined by 
other considerations. In particular, since lon is the small- 
est eigenvalue, the corresponding eigenfunction must be 
the translation mode of the defect that does not affect 
the energy. Therefore u>n = and the translation mode 
is known once we know w^: 



it 



a 

w N 



(16) 



where we have dropped the subscript oiia«. The profile 
function is: 



00 (x) = a [ — 
J w N 



(17) 



and Eq. (fTT7f> gives: 



V(<h) 



2w% 



(18) 



s(0o) 



The construction of U n is iterative and one must start 
with the highest bound state with eigenvalue To 
find w\ we set Uq = ojq to be a constant which will be 
determined by other considerations. Now we illustrate 
this scheme in a few cases. 



If the defect has only one bound state, it must be the 
translation mode. Since translations do not cost energy, 
the eigenvalue, w\, vanishes. So we need to solve Eq. fTJUl 
with u>i = 0: 



- w'l + UqWx = 
with Uq = uJq. After some manipulations 5 : 
ft = — ujq tanh(wox) 

and 

Ui(x) = f{ + fl = - 2sech 2 (^ x)] 

Now we can determine u>o by the requirement that at 
spatial infinity, the excitations are the particles in the 
trivial vacuum. If we denote the masses of these particles 
by to, this means that 

U(oo) = rn 2 



(19) 
(20) 
(21) 



Hence: loq — m and 

U(x) = to 2 [1 — 2sech 2 (ma;) 
The translation mode, tpt, satisfies 
- $ + U{x)i> t = 

and, in fact, the solution is 

ex. 



— = a scch(ma;) 



(22) 
(23) 
(24) 

(25) 



where a is a constant. Note that, in the present con- 
text, there is no requirement that i/jt be normalized as a 
wavefunction. 

Therefore the profile function is 



f , , / \ 2a 
• I ax sech(m.Tj = — tan 
J m 



tanh 



mx 



In other words, 



tanh 



mx 



— tan 



2a 



(26) 



(27) 



Since here the Bogomolnyi equations can be used, 
Eq. (|1C)|> with some algebra gives: 

„,2 



cos 



2mc 



1 



(28) 



This is the sine-Gordon potential. 

The potential still contains the unknown parameter a. 
To fix this parameter we could use some other property 
of the defect, for example, its total energy. 

This completes the global reconstruction of the poten- 
tial in the single bound state case. 



4 Note that the function w n (x) is not required to satisfy vanishing 
boundary conditions at infinity. In fact, by examining the dif- 
ferential equation 1131 in the asymptotic limit, it can be argued 
that ui n (±oo) will be divergent. 



We are closely following Ref . 7] . Note that the eigenvalues in 
are — k£ whereas we have taken them to be 
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B. Two bound states 

Now there are two eigenvalues ojf and wf . The trans- 
lation mode is always the lowest eigenvalue and hence 
u>2 = 0. So we first need to find a potential, U\(x) that 
contains the u)\ mode. 

Following the recipe given in the previous section, we 
have 



ui = fl + f{ + <4 

and we need to solve 



(29) 



(30) 



where, as before Uq = co 2 . This is exactly the single 
bound state problem that we solved in the previous sub- 
section if we replace w 2 with v 2 = lu{ 



and 



fx = —vtaxxhlyx) 



Ux = v 2 [l- 2sech 2 (i/x) 



jf. Therefore, 
(31) 



(32) 



Now we use the second eigenvalue. The potential U2 
is given by 

t/ 2 =/f +i2+^2-/f +/2 (33) 

where we have used u> 2 = since this mode must be the 
translation mode. So we need to solve: 

- w 2 + {v 2 [l - 2sech 2 (ra)] + uj 2 }w 2 = (34) 

We can rescale z = vx and bring the equation to the 
form: 



d 2 w 2 
dz 2 



+ {A + 2sech 2 (z)}u;2 = 



where: 



A= - 



1 + ^ 



(35) 



(36) 



Eq. has been solved in Ref. (see Section 6.3, 

page 768) and the solution is given in terms of hyperge- 
ometric functions: 



w 2 = m-m K/2 



2 Fx(K + 2,K-l;K + l;0 



where, 



and 



+ 2 Fx(K + 2,K-l;K+l;l-0 



_ 1 + tanh(z) 



K = -A 



(37) 



(38) 



(39) 



Note that we have fixed the ratio of the two linearly in- 
dependent solutions in Eq. (|37|l by imposing the require- 
ment that w 2 have even parity under z — > — z, which is 
the same as requiring W2 (£) = ^2(1 — ^). (The over- 
all normalization of W2 is unimportant.) Now Eq. i|18|) 
immediately gives the field theoretic potential though 
the expression is still implicit since the inverse function, 
x(<f>o), needs to be determined after doing the integral in 
Eq. fTfl . 

Let us now look at the special case when oj 2 /v 2 = 3 
(K = 2). Then the solution for w 2 is given by elementary 
functions: 



With this 



w 2 (x) = cosh 2 (r/x) 



r dx 

i>o(x) = olv \ — = atanh(ra) 
J w 2 



(40) 



(41) 



where a is a constant of normalization. (For convenience 
we have explicitly included a factor of v in the normal- 
ization.) Therefore: 



V^o)^sech 4 (,)^(l-§ 



(42) 



Hence we have recovered the A</> 4 model in this special 
case. 

It appears remarkable that the simplest examples of 
one and two bound states have led to two well-known 
field theories, namely the sine-Gordon and \<fi theories. 
A partial explanation is to be found in section [51 where 
we show that whenever U(x) arises from a field theoretic 
kink, it must necessarily be of the supersymmetric form 
(Eq. pi|l). So the recipe to reconstruct U(x) is ideally 
suited to the present problem. 



C. Three or more bound states; higher dimensions 

The reconstruction can be continued to the case of 
three of more bound states. However, the general re- 
construction becomes difficult to do analytically and nu- 
merical methods are desirable. Such methods have not 
yet been developed. 

The inverse scattering technique that we have used 
only works in one spatial dimension. Problems in more 
than one spatial dimension require other inverse scatter- 
ing techniques. However, if the problem is spherically 
symmetric, the problem effectively reduces to one spa- 
tial dimension and the reconstruction outlined above can 
be applied. (In Ref. Q the authors applied the inverse 
scattering technique to S-wave bound states of quarks.) 
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IV. APPLICATION TO COSMOLOGICAL 
INFLATION 

A. Inflation and defects 

If there was an extended but limited period of supcrlu- 
minal expansion in the early universe, light scalar fields 
would get excited by the spacetime expansion and create 
adiabatic perturbations on superhorizon scales. The pre- 
diction of adiabatic perturbations matches observations 
of anisotropies in the Cosmic Microwave Background Ra- 
diation (CMBR) and is the reason for the current confi- 
dence in the cosmic "inflationary" scenario. The mech- 
anism driving cosmic inflation is a scalar field called the 
"inflaton" . The dynamics of both the spacetime and the 
scalar field is dominated by the inflaton potential energy 
for the entire duration of cosmic inflation. Under cer- 
tain assumptions, this requires the scalar field to lie on 
top of a very flat potential. With time the scalar field 
slowly rolls along the flat potential, eventually reaching 
the steeper parts of the potential where its kinetic en- 
ergy becomes significant. Then inflation stops and the 
subsequent evolution of the scalar field leads to "reheat- 
ing" i.e. conversion of scalar field kinetic energy into 
ordinary matter in a thermal state 6 . 

One hopes that the scalar fields necessary for cosmic 
inflation will automatically arise in high energy particle 
physics models, such as Grand Unified Theories (GUTs) 
or string theory. On very general grounds, it is known 
that such theories contain magnetic monopoles and some- 
times other topological defects. If inflation and the de- 
fect are due to separate fields, the study of the defect 
core cannot yield information about inflation. For the 
global reconstruction scheme under discussion to work, 
the inflaton itself should play a role in the structure of 
the topological defect. 

One of the cosmological problems that inflation was 
designed to solve was the magnetic monopole overabun- 
dance problem that arose out of the marriage of GUTs 
and cosmology. If there is a period of inflation during or 
after the GUT phase transition, the magnetic monopole 
density would get diluted to insignificant amounts, leav- 
ing perhaps only a few in the entire visible universe. How- 
ever, subsequent work on inflation has shown that topo- 
logical defects can still be produced in significant num- 
bers toward the late stages |TTl HI HI 111 Il5j . If bound 
states on one of these defects can be studied, it would 
help in the global reconstruction of the inflaton poten- 
tial 7 . The other possibility of using the reconstruction 



scheme is if future particle physics experiments are able 
to produce magnetic monopole and antimonopole pairs 
(or closed walls or strings) in the laboratory 8 . The excita- 
tions of these monopoles could then be studied. There is 
also the possibility that duality holds in particle physics 
and the known particles (e.g. the electron) may be du- 
alized magnetic monopoles [Ty| . The internal excitations 
of the electron, would then be excitations of a magnetic 
monopole. Perhaps these excitations can tell us some- 
thing about inflation. 

Earlier work on reconstructing the inflaton potential 
(see e.g. [15) worked under the assumption that the 
CMBR anisotropies are generated by quantum fluctua- 
tions of the inflaton. This results in a local reconstruc- 
tion of the potential, the limited range of scales observed 
being explained entirely by a very small portion of the 
entire inflaton potential. Here wc will consider signa- 
tures of global features of the inflaton potential, such as 
flatness. 

Before we proceed to discuss signatures of inflaton po- 
tentials in the excitation spectrum, it is worth noting 
that inflation can occur within topological defects, so- 
called "topological inflation" We will, however, 
not consider this possibility since, if a defect inflates, it 
is not possible to find it or create it within our horizon 
while retaining predictability [2fJ. The condition for the 
topological defect to start inflating is that its width, 5, 
be larger than the horizon size, l/H, corresponding to 
the energy density within the defect. As we will show 
below (Eq. (|45|l l. in the case of domain walls with very 
flat potentials, 5 ~ r}/ y/2V(0), where rj is the change 
in cj) across the wall. The cosmological equations give 
H 2 = 8"7rGU(0)/3. Hence the condition for topological 
inflation is 7/ > rap where r\ is the vacuum expectation 
value of the field and rap is the Planck mass. Hence we 
will be restricted to r\ <C rap and we will ignore gravita- 
tional effects. 



B. Properties of the excitation spectrum 

The reconstruction of the inflaton potential can be car- 
ried out following the recipe given in Sec. [H] if the kink 
bound state spectrum is known. Here we will consider 
what is essentially the scattering problem for the infla- 
ton kink, namely the question: what signatures might we 
see in the bound state spectrum if the kink on hand is 
due to the inflaton? 

Consider a scalar field, cf>, with potential V((f>) that 



6 There are models of inflation that do not require a flat poten- 
tial. Instead they depend on a large initial value of the scalar 
field. Here we are limited to inflationary models that rely on the 
potential being very flat. 

7 Often the defects are due to topology present due to fields other 
than the inflaton. This significantly complicates the reconstruc- 
tion though, in principle, if there is variation of the inflaton 



within the defect, a reconstruction would still be possible. 
If inflation occurred at a relatively low energy scale - say some- 
what larger than the electroweak scale - this possibility is easier 
to envision and might even be realized at energies available at 
the Large Hadron Collider (LHC). However, the production of 
solitons by scattering particles is expected to be difficult because 
a large number of particles are simultaneously involved. 
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FIG. 1: Sketch of an inflaton potential. A crucial feature is 
the broad flat region from —771 to +771. The width of the 
curved region from ±771 to ±77 is assumed to be much smaller 
than 771. 



is invariant under the Z2 transformation <fi —* —cj> i.e. 
V(— 4>) = V(+(f>). Let the true vacua be given by (f> = ±77. 
Then there will be a domain wall solution across which 
<j) will change from —77 to +7/. We are interested in the 
case when <j> is also an inflaton with very flat V(<j>) for 
4> € (—77, +77). Here we will show that the flatness of the 
inflaton potential V((j>) has definite predictions for the 
spectrum of excitations of the domain wall 9 . 

A flat potential is drawn in Fig. ^ In the discussion 
below we will assume that (77 — 771) <C Tjx, so that the 
transition from flat to curved potential occurs relatively 
quickly 10 . A rough approximation to the potential in 
the interval (—77, +77) is given by a top hat with some 
modification near the curved regions in the vicinity of 
4> = ±771. The shape of the potential for </> > 77 will not be 
crucial for us but we will use the fact that the curvature at 
the global minimum is +to 2 where m is the mass of small 
excitations infinitely far away from the domain wall. The 
second derivative of the potential V"((f>) is sketched in 
Fig. [21 Note that the horizontal axis is <j> in this plot. 
The potential U(x) that determines the bound states is 
given by: 

U(x) = V"{cj> {x)) (43) 

where <fro(x) is the domain wall solution. (The problem 
is effectively 1+1 dimensional and hence we suppress de- 
pendencies on the y and z coordinates.) For the top hat 
potential: 

{—57 x < -8/2 
i]{2x/S) -6/2 <x<5/2 (44) 
+77 x > +6/2 



9 In non-Abclian field theories, domain wall solutions are much 
more complex. For ex ample, domain walls in SU(N) X Z2 have 
been discussed in Ref . |2 lL |22| . We will only discuss the simplest 
case of a Z2 kink with a single scalar field in this paper. 

If this assumption is not valid, the discussion below of the bound 
state spectrum will need to be modified. 



VOW 




FIG. 2: Sketch of the second derivative of the inflaton poten- 
tial. 

where 5 is the thickness of the defect and is determined 
by the Bogomolnyi equation © as: 

5 = V ^W) (45) 

Therefore the potential U(x) has the shape shown in 
Fig. |3 The broad flat bottom of the potential in the 
central region and the asymptotic behavior are generic to 
the inflationary models we are considering. The details 
of the transition regions near |x| ~ 6/2 may be model- 
dependent since they depend on the transition from flat 
V{4>) where the field rolls slowly to the curved part of 
V{4>) where reheating starts. The general features of 
U(x) are that is has a very broad central region where 
U(x) ~ 0, then a dip, then a rise that occurs in a width 
that is much smaller than S, and finally an asymptotic 
plateau of to 2 . 

We know that the lowest bound state is the translation 
mode and has vanishing energy. In fact, Eq. (|44|l can be 
used to find the translation mode tp t for the top hat case: 

_ dfo _ J 2 v /6 -6/2 <x<6/2 
m • ' dx \ otherwise 

Since the translation mode is the lowest energy bound 
state - an eigenstate with lower eigenvalue would signal 
an instability of the defect - all bound states have posi- 
tive energy and lie above the top of the double well struc- 
ture in the transition regions. So we expect the higher 
bound states to be relatively insensitive to the details 
of the transition region and mimicking U(x) by a finite 
square well potential may be a reasonable way to start a 
first analysis. 

The finite square well potential is analyzed in virtually 
every quantum mechanics textbook (e.g. |2^|). If Uoc 
denotes the depth of the square well and 8 its width, the 
number of bound states of a particle of mass fi is given 
by: 

N b ~ V^oc<5 2 (46) 

In our case (see Eq. J2J), 2/i = 1, Uoo — ™ 2 , and 8 is 
given by Eq. 14511 . Therefore: 

N b ~ -£2= (47) 
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U(x) 





m 2 






-8/2 


I ~~^\ 1 


8/2 x 



FIG. 3: Sketch of U(x) 



The three parameters m, r\ and V(0) occurring in this 
formula are independent: m 2 is the curvature in the true 
vacuum, V"{rf)\ rj is the position of the true vacuum; 
V(0) is the height of the flat part of the potential. If all 
these parameters were constrained by inflationary cos- 
mology, we would have some bounds on Nf,. However, 
there are no useful constraints on the parameter m and 
no model-independent bounds on can be derived. So, 
to get an idea of the range of Nf,, we work out its value 
in the case when V(<p) is of the Coleman- Weinberg form 
|24| with vanishing curvature at <f> = [25|: 



v{<t>) 



Wc then find: 



B 



In 



(48) 



(49) 



Hence we expect that an inflaton domain wall will have 
at least a few bound states. 

As is well-known [23], a special property of the infi- 
nite square well potential is that its eigenvalue spectrum 
is proportional to n 2 where n — 1,2,... is an integer 
that labels the eigenstates starting with the ground state. 
The same property holds for the low lying eigenstates of 
the finite square well 11 . If the potential is shaped like 
~ \x\ b , then it can be shown that E n oc n 2b ^ b+2 ' > (see 
Appendix Q. This dependence shows that observations 
of the spectrum of bound states on a domain wall can 
be used to find the power b and hence the shape of the 
potential. In the inflationary case, the width of the well 
is much larger than the distance over which the sides of 
the well get to their asymptotic levels. Therefore, the 
sides of the well are very steep (b ^> 1) in relation to the 
width of the wall and so the bound state spectrum w 2 
should be proportional to ~ n 2 . 

Qualitatively similar arguments may be used in the 
three dimensional case, in the case when the bound 
state spectrum on a monopole is known. If the po- 
tential, J7(x), is spherically symmetric, the density of 
energy eigenstates with fixed total angular momentum 



can once again lead to information about the flatness of 
the inflaton potential. For example, consider the s-wave 
states in a spherically symmetric potential U(r). The 
Schrodinger equation for the n th eigenstate radial wave 
function R n (r) = u n (r)/r reduces to: 



2/j, dr 2 



r + U(r)u n = E n0 u n 



(50) 



11 The exact eigenvalues depend on solutions of transcendental 
equations and must be found numerically. 



where the energy eigenvalue E n Q carries the I = label. 
As in the one dimensional case, here too we expect E n Q tx 
n 2 when U(r) is due to an inflaton. 

V. CONCLUSIONS 

We have described a "recipe" for recovering the poten- 
tial in a field theory, V((j>), starting with the bound state 
spectrum on a topological defect. An important aspect of 
the reconstruction discussed here is that it is "global" - 
the whole potential is reconstructed and not just a small 
part of it. As specific examples, we have applied the 
recipe to the case of one and two bound states on kinks. 
In the one bound state case, the recipe yields the sine- 
Gordon potential and in the two bound state case with a 
specific set of eigenvalues we obtained the X(f> 4 potential. 

While implementing the recipe we have relied on the in- 
verse scattering method based on supersymmetric quan- 
tum mechanics |7j. The method yields reflectionless po- 
tentials U(x). Other inverse scattering methods may also 
be used and, in general, they will lead to different U(x) 
with the same bound state spectra. The non-uniqueness 
of the reconstruction may be reduced by further inclusion 
of scattering data and perhaps by using other system- 
specific information. 

Our recipe for the reconstruction of V(cj>) is sufficiently 
involved that, except in the simplest situations, it will 
have to be implemented numerically. For example, even 
in the two bound state case, the solutions of the differen- 
tial equations in the inverse scattering problem are given 
in terms of hypergeometric functions, making it very hard 
to analytically implement the recipe. The reconstruction 
in 3+1 dimensions is technically even more challenging. 
However, we can expect the analysis to be close to the 
1+1 dimensional case when the problem is spherically 
symmetric. 

The results of this paper are relevant to any system 
in which topological defects occur. Hence we can ex- 
pect applications to condensed matter systems, parti- 
cle physics, and cosmology. It would be interesting to 
try out the reconstruction in condensed matter systems 
where topological defects are readily available. (Or in nu- 
clear physics to the extent that nuclei can be modeled by 
Skyrmions [2ti 1271 l2q .1 The reconstruction would yield 
a Landau-Ginzburg type of effective interaction potential 
but would not yield, at least directly, information about 
the microphysical interactions between the fermions. The 
application to particle physics and cosmology is futuris- 
tic since topological defects are theoretically expected in 
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these settings but have not yet been experimentally dis- 
covered or observed. Just as in the condensed matter 
case, the scalar field need not be fundamental even in 
the particle physics context. 

A novel application of the reconstruction recipe is in 
the context of inflationary cosmology. If the inflaton 
vacuum manifold has suitable topology topological de- 
fects in the inflaton field will exist. The spectrum of 
bound states in these defects will reflect the properties of 
the inflationary potential. Reversing the argument, since 
we know that the inflationary potential must have cer- 
tain properties, the bound state spectrum must also have 
some characteristics. We have discussed these character- 
istics as a way to probe the inflationary scenario, which 
is hard to do otherwise • If future investigations dis- 
cover a scalar field with an extremely flat potential, that 
scalar field will be a prime suspect to be the inflaton, and 
the defect with its characteristic bound state spectrum 
will be a "smoking gun" from the shot that was fired ten 
billion years ago. 

If the topology of the inflaton vacuum manifold is triv- 
ial, no defects will exist and the reconstruction scheme 



discussed here will not be useful for inflationary cosmol- 
ogy. If, however, inflaton topological defects do exist, 
experiments may become feasible in the future that can 
directly probe cosmology in the laboratory. 
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APPENDIX A: CHECK OF THE INVERSE 
SCATTERING EQUATIONS 

Here we give a check of the iterative scheme for the 
inverse scattering method described in Ref. 0- 

Suppose we are given a potential U n -\(x) that contains 
n — 1 of the highest bound states. Then we construct the 
functions f n (x) from the equation: 

fn - fn + Un-1 = l«£ (Al) 



9 



This non-linear equation is simplified by setting 

w' 

fn = (A2) 

using which we get the linear equation 

- < + U n -iw„ = Lo 2 n (A3) 
Finally we construct 

U n {x)=f n + fl+ul (A4) 

The claim is that U n has an eigenstate with eigenvalue 
w 2 in addition to all the other higher eigenstates. 

That U n admits a state with eigenvalue w 2 , can be 
shown explicitly. The state is given by: 

fl>n = — (A5) 

Then it easy to check that 

- 1p'£ + U n 1p n = UJn^n (A6) 

So, indeed, the potential U n has an (explicitly con- 
structed) eigenstate with eigenvalue to n . 

For the iterative procedure to work, we also need to 
show that U n admits eigenstates with the higher eigen- 
values {w 2 } for i = 1, ...,ti — 1. This is seen as follows. 
U n and U n -i are "partner" potentials. In other words, 
we can write the n th Hamiltonian as: 

H = {-dl + U n ) = A+A + uj 2 n (A7) 

where 

A = —d x + f n , A+ = +d x + f n (A8) 

Then, the partner Hamiltonian is: 

H_ = (-dl + Un-i) = AA+ + u? n (A9) 

It is now easy to show that if tp satisfies 

Hip = Eip (A10) 

then Aip is an eigenstate of H- with the same eigenvalue: 

H_(Aip) = E(Aip) (All) 

Thus, H and H- have a common eigenspectrum, except 
for the lowest state of H that satisfies Aip = 0. Hence 
U n has states with eigenvalues tof, i — 1, ...,n — 1 since 
these are also the eigenvalues of states in U n ~i, and then 
it has one extra eigenstate and this has eigenvalue uj n . 

APPENDIX B: FORM OF DOMAIN WALL 
FLUCTUATION POTENTIALS 

Here we show all that all kink potentials (U) have 
the supersymmetric form. We know that every kink has 



U(x) 



_ -w m 2 


+w 
















/ E " 


-8/2 


\ — ^\ / 


8/2 x 



FIG. 4: Sketch of U(x) together with energy level E n . The 
shaded region is the interior of the square well potential, U sq , 
that is used to find the approximate dependence of E n on n. 

a translation mode which does not change the energy. 
Therefore, the Schrodinger equation (Eq. (J2Jl) gives 

U( X ) = f (Bl) 

where ipt(x) is the translation mode. This form of U(x) 
can be rewritten as: 

U(x) = f' + f 2 (B2) 

where / = [ln(^t)]', and hence U is of the supersymmet- 
ric form (Eq. ((TTj) ) |30|. 

This argument can also be extended to higher dimen- 
sions with derivatives replaced by higher dimensional 
derivatives. Then we find: 

U = V • f + f 2 (B3) 

where f(x) = V(ln-0t) is a vector valued function and 
i// t (x is a translation mode. The Hamiltonian is: 

H = A+ - A (B4) 

where 

A=-a + f (B5) 

If tp is a multi-component scalar field, or a collection 
of several scalar and gauge field fluctuations, in several 
dimensions, there are circumstances in which U is still 
supersymmetric. Label the many components of tp by the 
index i. So tp may be thought of as a column vector with 
components tpi- Now in higher dimensions, there will be 
several zero modes. For example, translations along any 
dimension will be a zero mode. So label the zero modes 
by the index a and denote the translation modes by £ a . 
Then consider the matrix M whose components are 
We can show that U is supersymmetric if we assume that 
M is a square matrix that is invertible. This follows 
because now: 

[-V 2 + U(x)]M = (B6) 

(U itself is a matrix potential.) Then it is easy to show 
that 

U(x) = V • F + F 2 (B7) 

where 

F = (VM)lVr 1 (B8) 
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APPENDIX C: SPECTRAL PROPERTIES OF 
BOUND STATES IN 1+1 DIMENSIONS 

We now find the connection between the shape of the 
potential and the dependence of the n th eigenvalue, E n , 
on n for the specific class of potentials: 



U(x) = a\x\ 



(CI) 



where a and b are parameters. 

We use the WKB approximation {2^. The quanti- 
zation condition for a particle of mass fj, — 1/2 in the 
potential U{x) is 



+w(E n ) 



■j(E n ) 



dx\/ E n — U (x) = n 



(C2) 



where E n stands for w 2 in the notation of Sec.[n] Insert- 
ing the class of potentials in Eq. (|C1|I . we find: 



(C3) 



Next we conside the WKB method in the context of 
the potential U(x) redrawn here in Fig.^J 
Now, 

p+w{E n ) 

n = I dx\/E n — U(x) 

J-w{E n ) 
r+w(E n ) 



-w(E n ) 



dx^(E n -U sq )-{U -U sq ) (C4) 



U sq denotes the infinite square well potential of width 
2w(E n ) (see Fig. EJ. Our assumption is that the domi- 
nant contribution to the integral comes from the E n — U sq 
term and that the U(x) — U sq term can be ignored. Then 
the integral is trivial to do, resulting in n — \/E n 2w(E n ). 
Now 2w{E n ) - 5. Therefore E n ~ n 2 . 

If E n is very small, our assumption that U (x) — U sq 
can be ignored will not hold. For reasonably large E n 
we can expect it to hold, though E n should still be much 
smaller than m 2 - the asymptotic value of U - so that 
the finite depth of the well does not play a role. 



